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The strange quark scalar content plays an important role in both the description of nucleon struc- 
ture and in the determination of dark matter direct detection cross sections. As a measure of the 
strange-quark contribution to the nucleon mass, the strange-quark sigma term (<7 S ) provides im- 
portant insight into the nature of mass generation in QCD. The phenomenological determination 
of a s exhibits a wide range of variation, with values suggesting that the strange quark contributes 
anywhere between and more than 30% of the nucleon mass. In the context of dark matter 
searches, coupled with relatively large Higgs coupling to strangeness, this variation dominates 
the uncertainty in predicted cross sections for a large class of dark matter models. Here we report 
on the recent results in lattice QCD, which are now giving a far more precise determination of 
a s than can be inferred from phenomenology. As a consequence, the lattice determinations of 
a s can now dramatically reduce the uncertainty in dark matter cross sections associated with the 
hadronic matrix elements. 



The 30 International Symposium on Lattice Field Theory - Lattice 2012, 
June 24-29, 2012 
Cairns, Australia 



* Speaker. 



© Copyright owned by the author(s) under the terms of the Creative Commons Attribution-NonCommercial-ShareAlike Licence. 



http://pos.sissa.it/ 



Strange quark content of the nucleon and dark matter searches 



R. D. Young 



1. Introduction 

The modern picture of the universe suggests that the ordinary matter component in the energy 
composition of the universe is only about 4% — with the remainder comprised of some form 
of cold, "dark matter", which clusters around ordinary matter and an even larger "dark energy" 
component. While little is known about the physical nature of dark energy, there is strong evidence 
that suggests we are nearing the discovery phase in the identification of dark matter. Supposing 
dark matter to take a particle-like form, the relevant mass scale for such particles are most likely to 
be within reach of the LHC. 

For a general spin-independent interaction of a WIMP with a nucleus, the low-energy limit 
reduces to a scalar contact interaction, and hence sensitive to the qq matrix elements within a nu- 
cleon. And for Higgs-dominated exchange, these couplings are proportional to the corresponding 
quark mass and are hence sensitive to the nucleon sigma terms. The sigma terms therefore become 
of principal uncertainty in the predicting the cross sections associated with any candidate model 
of dark matter. The significance of the uncertainty in these hadronic matrix elements has been 
highlighted in a range of dark matter models, see for example Refs. [1, 2, 3, 4, 5]. 

Before going on to discuss the sigma terms, and the strangeness scalar content, it is worth 
recapping some recent advances in resolving the strange quark vector form factors in the nucleon. 
Early estimates of the strangeness electromagnetic currents had suggested that they could be rela- 
tively large [6]. After a dedicated experimental effort, parity-violating electron scattering measure- 
ments [7,8,9] have revealed that the strange quarks contribute much less than originally suggested. 
This finding is also supported by lattice-based phenomenological estimates [10, 11] and recent di- 
rect lattice QCD simulation results [12, 25]. For a recent review of these latest revelations, see 
Ref. [13]. 

After setting some general notation, the phenomenological determination of the sigma terms 
are reviewed in Section 2; a summary of the latest lattice QCD results are reported in Section 3; 
the significance of these results in the context of dark matter searches are discussed in Section 4; 
follwed with a summary in Section 5. 

1.1 Notation 

Of primary interest here are the scalar nucleon matrix elements, where we'll use the notation 
for the light- and strange-quark matrix elements 

C/ = mi{N\uu + dd\N) , a s = m s (N\ss\N) , (1.1) 

with the average light-quark mass given by mi = (m u -\-m c i) /2. 

A commonly reported measure of the strangeness sigma term is through the y-parameter, 

y _ 2(N\ss\N) = 2m, a s 
y ~ (N\uu + dd\N) m s o, ' 

2. Phenomenological Determination 

By the very nature of the weak coupling of the Higgs to the low-energy sector of the Standard 
Model, the sigma terms are essentially impossible to measure directly. Fortunately, G/ can be in- 
ferred through a chiral low-energy relation, where the amount of explicit chiral symmetry breaking 
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can be related to pion-nucleon scattering. In particular, the light-quark sigma term can be extracted 
through measurement of the Born-subtracted, isoscalar amplitude L n ff(t). 

Both Ljzn and a vanish with the quark mass, but importantly they become equal as the chiral 
limit is approached. That is, the leading dependence of the quark mass is the same, with the leading 

3 /2 

difference being &\m, ). The limited knowledge of this difference term can further be reduced 
by moving to the unphysical kinematic point t = 2m\ (the Cheng-Dashen point [14]), where the 
remainder is ff(mf) [15], being defined by 

A R = Z m (t = 2ml) ~ °l{t = 2»4) = 0{m 2 i) . (2.1) 

Here, the scalar matrix element has been extended to non-zero momentum transfer, with the usual 
sigma term corresponding to the t — > limit, C/ = Gj(t = 0). An early calculation of the remainder 
term has determined an estimate A# ~ 0.35 MeV [16], later followed by an updated value A# ~ 
2MeV [17]. 

Extraction of the sigma term at t = then requires the determination of the form factor cor- 
rection, defined by 

A a = ai(t = 2ml)-ai(t = 0). (2.2) 

Through dispersion relations, this form factor correction has been estimated to be A CT ~ 15 MeV 
[18]. 

With the required theoretical corrections under reasonable control, the pion-nucleon sigma 
term can then be extracted from the experimentally determined as extrapolated to the unphys- 
ical Cheng-Dashen point. In summary, <j/ is determined by 

a, = Z m (t = 2ml) ~ A R ~ A a- (2-3) 

Following this outlined technique, analysis of 7tN scattering data gives the Gasser-Leutwyler- 
Sainio (GLS) value C/ = 45 ± 8 MeV [16], or a somewhat larger value from the George Washington 
University/TRIUMF (GWU) group analysis o } = 64 ±7 MeV [19]. 

In addition to these benchmark calculations of GLS and GWU, an analysis based on a covariant 
baryon chiral perturbation theory has recently been reported by Alarcon, Martin-Camalich & Oiler 
(AMO) [20]. Here the low-energy 7iN scattering phase shifts are fit directly within the effective 
field theory framework. The analysis determines the relevant low-energy constants necessary for 
the extraction of the sigma term by the Hellmann-Feynman theorem (discussed below). Together 
with estimates for the higher-order terms, the sigma term is reported to be a/ = 59 ± 7 MeV, lying 
between the two values reported above. 

Estimating the strange-quark sigma term is significantly more challenging. In this case, the 
strange quark is too heavy to reliably truncate the low-energy relation at low order. As an alter- 
native, the conventional approach has been to estimate a s by studying the SU(3) breaking among 
the baryon octet [21, 22]. Here, the baryon mass splittings can be used to constrain the non-singlet 
combination 

a = mi{N\uu + dd-2ss\N) . (2.4) 
To leading-order in the quark masses, Oq can be estimated from the physically observed spectrum 

by 

Ob ^ — (Mz + M^- 2M N ) ~ 24 MeV . (2.5) 

m s — mi 
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Figure 1: The determination of o s , from the best-estimate of at) = 36 ±7 MeV [22], has a high degree of 
sensitivity to the determination of pion-nucleon sigma term. The three values for a\ are as discussed in the 
text. 

By incorporating the higher-order terms in the quark mass expansion, Borasoy and MeiBner have 
determined an improved estimate Ob = 36 ±7 MeV [22]. With an estimate for On. the strangeness 
sigma term is then given by 



Being multiplied by the large quark mass ratio, this method leads to a value for o s that is acutely 
sensitive to the difference between C\ and On- Further, given the limited precision available to On. 
even a perfect determination of 07 leaves a residual uncertainty in a s of order 90 MeV. For the three 
determinations of C\ discussed above, Figure 1 displays the broad range of possible o s values. 

Given the difficulty in extracting a precise determination of a s from phenomenology, there is 
significant scope for lattice QCD to provide a meaningful constraint on this nucleon matrix element. 
In addition, there is also the potential for lattice simulations to shed light on the phenomenological 
determination of a). 

3. Lattice QCD 

Within the framework of lattice QCD, there are two main methods used in the extraction of 
the sigma terms. These divide into the explicit evaluation of the scalar qq matrix element by 3- 
point function methods or by invoking the Hellmann-Feynman relation through the study of the 
quark-mass dependence of the nucleon mass. 

In the first technique, appropriate ratios of 3-point and 2-point correlation functions are formed 
to isolate the matrix elements of interest. The evaluation of the correlation functions involves two 
distinct forms of Wick contraction, as depicted in Figure 2. In particular, there are contributions 




(a, - oo) • 



(2.6) 
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1=1 B 

Figure 2: Schematic diagram displaying the two topologically distinct contractions in the evaluation of 
the scalar 3-point function. The standard jargon is that the left is a connected insertion and the right a 
disconnected insertion. 



involving quark line connected or disconnected operator insertions. While standard methods typ- 
ically yield strong signals for the connected correlation functions, it is well-known that the dis- 
connected insertions are notoriously challenging, see eg. Refs. [23, 24, 25]. This of particular 
significance for the strange-quark matrix elements in the nucleon, which are purely disconnected. 

A common alternative to the three-point method, is to determine the qq matrix elements by 
differentiation with respect to the quark masses, where the Hellmann-Feynman relation gives [26, 
27, 28] 

dM N 

O q =m q - . (3.1) 

am q 

Here the Gell-Mann-Oakes-Renner (GOR) relation [29] is commonly imposed such that Mjy is 
expressed as a function of the squares of meson masses. One way to do this is to write Mn = 
M^(m^,m|-), where m\ = m 2 K — m\/2 is the projection of the square of the kaon mass onto the 
SU(2) chiral limit (m/ — > 0). With such a formulation, the sigma terms (for 2+1 -flavour simulations) 
are easily written as 

2 dM N „ 2 dM N 
Oi = m % 3— J- , a v = m K —j- . (3.2) 
am l n dm A K 

The main challenge of this approach is the difficulty in reliably parameterising the quark-mass 
dependence over a range of light and strange quark masses. In particular, it is only in recent years 
that there have been large scale numerical simulations of baryons with 2+1-flavours of dynamical 
quarks, eg. [30, 31, 32, 33, 34, 35, 36]. In constraining the two-dimensional parameter space, it 
is also the case that typical lattice trajectories in the m\-m s plane approach the physical point for 
approximately constant m s ; though the QCDSF-UKQCD are currently pursuing an approach which 
keeps the singlet combination (2m/ +m s ) a constant [36]. 

A summary of the progression of lattice results is displayed in Figures 3 and 4. With differing 
degrees of analysis into the various uncertainties of the calculations, it would be difficult to formu- 
late any rigourous aggregates. Nevertheless, two clear features are emergent. For G/, the values 
revealed in lattice simulations are compatible with the range of phenomenologically determined 
values. Secondly, the modern values of G s are all at the lower end of the possible values suggested 
by Figure 1. To highlight the current status, Figure 5 shows a close-up of the recent determinations 
of <7... 
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Figure 3: Light-quark sigma term results based on lattice QCD. The colours denoted the number of dy- 
namical flavours of quarks: green is Nf = 0, blue Nf = 2 and red Nf > 2+ 1. References: Fukugita et 
al. [37], Dong et al. [38], SESAM [39], Leinweber et al. (2000) [40], Leinweber et al. (2003) [41], Procura 
et al. (2003) [42], Procura et al. (2006) [43], ETM [44], JLQCD [45], QCDSF (201 1) [46], QCDSF (2012) 
[47], Young & Thomas [48], PACS-CS [49], Martin-Camalich et al. [50], Diirr et al. [51], QCDSF-UKQCD 
[52], Shanahan et al. [53], Ren et al. [54]. 

4. Dark Matter 

The smaller values of a s revealed in the recent lattice studies are particularly significant in 
the context of the direct search for dark matter. The most precise limits on WIMP-nucleon cross 
sections are being constrained by the XENON 100 Collaboration, with the latest update placing an 
upper bound on the cross section of less than 10~ 44 cm 2 over a wide range of WIMP masses [61]. 
Figure 3 of [61] suggests these limits are continuing to reduce the parameter space of potential 
supersymmetric candidates for dark matter. 

The XENON 100 Collaboration results are plotted against predicted cross sections for some 
favoured supersymmetric models [62, 63, 64]. The predicted cross section rates are based on a 
determination of the strange quark sigma term, a s } as outlined in Section 2. Hence a s in these 
studies exhibits the extreme sensitivity to C\ displayed in Figure 1 . 

As the WIMP-nucleon interactions are largely Higgs-coupling driven, the difference between 
a small and large o s can have a dramatic influence on the predicted cross sections. This is high- 

Or in an alternative common notation, fj s = a s /M p , for the proton mass M p . 
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Figure 4: Strange-quark sigma term results based on lattice QCD. Colours as in Figure 3. Fukugita et 
al. [37], Dong et al. [38], Lewis et al. [55], SESAM [39], JLQCD (2008) [45], JLQCD 2010 [56], QCDSF 
[46], Young & Thomas [48], Toussaint & Freeman [57], Martin-Camalich et al. [50], Diirr et al. [51], 
QCDSF-UKQCD [52], Freeman & Toussaint [58], Shanahan etal [53], JLQCD (2012) [59], Ren etal [54], 
Engelhardt [60]. 



lighted in Figure 6, which shows how the predicted cross section for a particular constrained min- 
imal supersymmetric standard model (CMSSM) model 2 depends strongly on E^jy (with a s con- 
strained by the phenomenological On) [65] . In contrast, the displayed ellipse shows the range of 
predicted cross sections within the 95% confidence level interval of the lattice QCD determinations 
of a) and a s from Refs. [48, 57]. It should be stressed that the reduced variation in the cross section 
is a consequence of the increased precision in a s from lattice QCD input — which is not reliant on 
the propogation of the phenomenological uncertainty in Oq. 

Generic dark matter cross section packages, such as micrOMEGAs [68], have been designed 
to take as inputs 0/ and On- With the improvement in lattice QCD results discussed above, it would 
be advantageous to see these packages reformulated to take oj and a s as inputs 3 . In the meantime, 
with cross section predictions based on do as an input, the reduction in uncertainty in a s may 
be equivalently stated as a reduction in a/ — Ob, cf. Eq. (2.6). A crude, yet conservative view of 



2 The figure displays the predicted cross-section for "model C", as one of a class of benchmark models proposed 
pre-LHC [66, 67]. 

3 Of course these are precisely the same thing with an appropriately included correlation coefficient. 
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Figure 5: Zoomed in graphic of Fig. 2 showing more recent results on a s . 
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Figure 6: For a pre-LHC CMSSM dark matter model, the predicted spin-independent cross section (<7si) 
shows a strong dependence on T. k n. This variation is a consequence of the large variation of O s , as con- 
strained by the phenomenological Go, as seen in Figure 1. The ellipse show the range of Osi using lattice 
inputs for L k n and a s [65]. 
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Figure 5, may suggest a value 4 

a,-a = — (40±30MeV) ~ 2.9±2.2MeV. (4.1) 

m s 

Already at this scale of precision, there should be a substantial reduction in the uncertainties in 
Csi associated with the hadronic matrix elements. Importantly, with the reduction in the hadronic 
uncertainty, any discovery of dark matter will have significantly more discrimination power among 
candidate models. 

5. Summary 

To summarise, an accurate determination of the strange quark sigma term is of principle im- 
portance in the reduction of hadronic uncertainties in the predicted dark matter cross sections for a 
wide range of models. In the determination of the relevant nucleon scalar matrix elements, lattice 
QCD simulations have made significant progress in recent years — particularly with the emergence 
of dynamical simulations with 2+1 flavours of dynamical quarks. As a confirmation of lattice meth- 
ods, it is reassuring to observe that the pion-nucleon sigma term appears to be compatible with the 
reliably determined phenomenological extraction. Further, recent determinations of the strange 
quark sigma term are significantly smaller than had previously been suggested — and lattice cal- 
culations are now at a far greater precision. With the prospect of a discovery of dark matter in 
the near future, it will be essential for lattice QCD simulations to further reduce these hadronic 
uncertainties. 
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